Electron cooling by diffusive normal metal - superconductor tunnel junctions 
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We investigate heat and charge transport in NN'IS tunnel junctions in the diffusive hmit. Here N and S are 
massive normal and superconducting electrodes (reservoirs), N' is a normal metal strip, and I is an insulator. 
The flow of electric current in such structures at subgap bias is accompanied by heat transfer from the normal 
metal into the superconductor, which enables refrigeration of electrons in the normal metal. We show that the 
two-particle current due to Andreev reflection generates Joule heating, which is deposited in the N electrode 
and dominates over the single-particle cooling at low enough temperatures. This results in the existence of a 
limiting temperature for refrigeration. We consider different geometries of the contact: one-dimensional and 
planar, which is commonly used in the experiments. We also discuss the applicability of our results to a double- 
barrier SINIS microcooler. 

PACS numbers: 74.45.+C, 74.50.+r, 74.40.Gh, 74.25.fc 
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I. INTRODUCTION 

The flow of electric current in NIS (Normal metal - Insula- 
tor - Superconductor) tunnel junctions is accompanied by heat 
transfer from the normal metal into the superconductorii^ 
This phenomenon arises due to selective tunneling of high- 
energy quasiparticles out of the normal metal which is induced 
by the superconducting energy gap. It is similar to the Peltier 
effect in metal-semiconductor contacts'* and enables refriger- 
ation of electrons in the normal metal. The heat current out 
of the normal metal (also refeiTed to as "cooling power") is 
maximal at a voltage bias just below the energy gap, eV < A. 
For eV > A both the current / through the junction and the 
Joule heating power IV strongly increase, rendering the cool- 
ing power negative. 

A micrometer-sized refrigerator, based on a NIS tunnel 
junction, has been first fabricated by Nahum et al. ^ The au- 
thors used a single NIS junction in order to cool a small nor- 
mal metal strip. Later Leivo et al^ noticed that the cooling 
power is an even function of the applied voltage, and fabri- 
cated a refrigerator with two NIS junctions arranged in a sym- 
metric series configuration (SINIS). This results in reduction 
of the electron temperature from 300 mK to about 100 mK, 
offering perspectives for the use of NIS junctions for on-chip 
cooling of nano-sized systems, like high-sensitive detectors 
and quantum devices.^ To enhance the performance of NIS 
microcoolers, it is important to understand possible limita- 
tions of the NIS refrigeration. 

Serious limitations of the cooling effect arise from the fact 
that nonequilibrium quasiparticles injected into the supercon- 
ducting electrode accumulate near the tunnel interface. ^-^ The 
consequences are the backtunneling of hot quasiparticles to 
the normal metal**-^, the emission of phonons (by the recom- 
bination of nonequilibrium quasiparticles into Cooper pairs) 
that partially penetrate the normal metal,-^ and the overheat- 
ing of the superconducting electrode.^ All these effects re- 
duce the efficiency of NIS refrigerators. This problem can be 
solved by imposing a local thermal equilibrium in the super- 
conducting electrode.^ So called quasiparticle traps /"-*' made 



of an additional normal metal layer covering the supercon- 
ducting electrode, remove hot quasiparticles from the super- 
conductor and are thus beneficial in this respect. 

However, there is a fundamental limitation for NIS micro- 
coolers. It arises from the intrinsic multiparticle nature of cur- 
rent transport in NIS junctions which is governed not only by 
single-particle tunneling but also by two-particle (Andreev) 
tunneling. The single-particle cuiTent and the associated heat 
cuiTent are due to quasiparticles with energies £ > A (com- 
pared to the Fermi level). At very low temperatures, single- 
particle processes are exponentially suppressed in the subgap 
voltage region eV < A, and the charge is mainly transferred 
by means of Andreev reflection of quasiparticles with ener- 
gies E < A. *^ *^ The Andreev current Ia does not transfer heat 
through the NS interface but rather generates the Joule heat- 
ing IaV which is deposited in the normal metal electrode — 
and dominates single-particle cooling at low enough tempera- 
tures. Thus the interplay between the single-particle tunneling 
and Andreev reflection sets a limiting temperature for the re- 
frigeration. 

The role of the Andreev cuiTent in the electron refrig- 
eration has been first theoretically analyzed by Bardas and 
Averin for the simplest model of the NIS microcooler - a one- 
dimensional constriction between the N and S reservoirs^^ 
assuming the constriction length to be much shorter than the 
coherence length. In experiment, the importance of Andreev 
processes in NIS microcoolers was first demonstrated by Ra- 
jauria et al,^'^ by using the theoretical estimations of the An- 
dreev cuiTent,*^ obtained within the tunnel Hamiltonian tech- 
nique, for interpretation of the experimental data. In this paper 
we present a quantitative analysis of heat transport in diffu- 
sive NIS tunnel junctions based on the solution of microscopic 
equations of nonequilibrium superconductivity. We consider 
the general case of arbitrary length of the normal wire, as 
well as of different possible geometries of the junction: one- 
dimensional (ID) junctions and planar junctions with overlap- 
ping thin-film electrodes, commonly used in experiments.'- 
We also discuss the applicability of our results to a double- 
barrier SINIS microcooler. 
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FIG. 1: One-dimensional (a) and planar (b) models of the NN'IS 
junction. The insulating barrier is shown by thick black line. 

The paper is organized as follows. In the next Section, we 
develop a theory for ID junctions. We start with a discus- 
sion of basic equations and adopted approximations, calculate 
the spectral characteristics of the junction using Usadel equa- 
tions, and finally obtain both the electric and the heat currents 
through the junction. In Sec.|lll]we extend this theory to the 
case of planar junctions. We discuss the results in Sec.lIVIand 
then consider possible extension of our theory to the case of a 
double-barrier SINIS junction in Sec. |V] Finally, we summa- 
rize the results in Sec.lVIl 

II. ID NN' IS JUNCTION MODEL 
A. Basic equations 

The model of the one-dimensional NN'IS junction is de- 
picted in Fig.lTJa). It consists of a voltage-biased normal me- 
tal reservoir (N) and a normal metal wire (N') of length L con- 
nected to a superconducting reservoir (S) through an insulator 
layer (I). We assume the NN' interface to be fully transparent. 

In our theoretical analysis, we consider the diffusive limit, 
in which the superconducting coherence length is given by 
expression ^'D/2A, where 2) is the diffusion coefficient 
of the normal metal (we assume kg — 1) and the elastic 
scattering length I ^ ^q. In this case, calculation of the elec- 
tric and heat currents requires solution of the one-dimensional 
Keldysh-Usadel equations'^ (see also the review'^) for the 
4x4 matrix Keldysh-Green function G{x,E) in the N' lead, 

[a-E , G] = iDdJ, J = Gd^G, ( 1 ) 

G=({ J^), G'^^ff-ff. (2) 

Here C5V is the Pauli matrix, dx = d /dx, g^-^ are the 2x2 
Nambu matrix retarded and advanced Green functions, and 
/ = /+ + CT^f- is the matrix distribution function (we use 
'check' for 4x4 and 'hat' for 2x2 matrices). In Eqs. ([T]l 
we neglect the inelastic collision term, assuming the length L 
of the N' lead to be smaller than the inelastic relaxation length. 

Equations ([T]i can be decomposed into the diffusion equa- 
tions for the Green functions, 

[a,E,g] = 1^9,1 J = gd,g, f = (3) 



and the equation for the Keldysh component G^, 

[(T,£, G^] = i'DdxF, F = fdxG'' + G^a,/. (4) 

Taking into account the normalization condition = 1, we 
parameterize the Green function by the complex spectral an- 
gle 0, 

g{x,E) = CT, cosh + /(7y sinh . (5) 

The electric and energy currents are related to the Keldysh 
component of the matrix current J as'^''§ii2i22 

I=— f I^dE, f EI+dE, (6) 

e Jo Jq 

/_ = iTr(J,/^=D_a,/_, I+=^TrF ^D+dJ+, (7) 

where g^ is the normal conductance of the N' lead per unit 
length, and D± are dimensionless diffusion coefficients, 

D = (l/4)Tr(l-(7,/o;|'^) =cosh2(Re0), (8a) 
D+ = (l/4)Tr(l~|V) =cos2(Im0). (8b) 

Here we expressed the advanced Green functions through the 
retarded ones using the general relation g^ = —O^.g^^'^^.r^ 
then omitted the superscript R. The quantity /+ has the mean- 
ing of the spectral density of the net probability current of 
electrons and holes, while /_ represents the spectral density of 
the electron-hole imbalance current responsible for the charge 
transfer (see the discussion in Refs. 21 and 18). 

Calculation of the electric and energy currents in Eqs. ^ 
involves two steps: first one has to solve the diffusion equa- 
tions (O for the spectral angle 0, and then to solve the kinetic 
equations (|4| to find the distribution functions. 

The expression for the heat current out of the normal metal 
reservoir (cooling power) in a diffusive NIS structure was sug- 
gested by Bardas and Averin in Re f. ilSl In contrast to the 
ballistic case (also discussed in Ref. [Tsl) . its spectral density 
contains several additional terms which are odd in energy and 
therefore vanish upon integration over energy. We propose 
another method for the derivation of the cooling power, which 
seems to be physically clearer and does not involve the non- 
physical, odd-in-energy terms. We define the heat generation 
in the reservoir through the work done by the applied voltage 
on the nonequilibrium quasiparticles coming to this reservoir, 
i.e., through the change of the kinetic energy Ej^ of the quasi- 
particles. We accept the definition E^ — E — e(p{x) for the 
electron-like and E'^=E + e(p[x) for the hole-like quasiparti- 
cles, where (p{x) is the electric potential (note that the quantity 
E is the total quasiparticle energy which is conserved during 
passage across the junction, in contrast to E/^). Along this line 
of reasoning, the heat generation in a given reservoir can be 
defined as the kinetic energy flow to that reservoir, 

4 (x) = ±NfS £° [E't ix)I'' + e'I (.^ )/''] dE. (9) 

We take x = and the minus sign for the left S reservoir, and 
x = L and the plus sign for the right N reservoir; Nf is the 
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electron density of states per spin in the normal state, S is the 
cross-sectional area of the junction, and the quantities I" and 
are the electron and hole probability flow densities, respec- 



tively. The expressions for F and were found in Ref.ullbv 
introducing the following parametrization of the matrix distri- 
bution function (see also Ref. 



/=l-2 
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The functions and n'' have the meaning of the electron 
and hole population numbers, respectively, and approach the 
Fermi distribution in the reservoirs. Then the electron and 
hole probability currents are defined as^ 

= (l/2)D(/+±/_) (10) 
= -{l/2)V[{D+± D^)d,n' + (D+ T D^)d,n^] . 

In the N reservoir (0 = 0, D± = 1), the currents /''■'' are 
naturally related to the electron and hole diffusion flows, 
= -Da,«f ''. Within the N' lead each current I'-^ gen- 
erally consists of a combination of both electron and hole dif- 
fusion flows, which reflects the coherent mixing of electron 
and hole states in the proximity region. Upon substitution of 
Eq. ( [Tol l into Eq. (|9]l, using the relation = le^Np'OS, we 
obtain the well-known equation for the heat current out of the 
normal metal reservoir (cooling power). 



-4(L) = -/y-e. 



(11) 



In the case of a NIN (Normal metal - Insulator - Normal metal) 
structure, the heat generation in both reservoirs was found to 
be equal to IV /2^¥ox the NIS structure, it is the imbalance 
between the kinetic energy flows to the N and S reservoirs 
that leads to the cooling effect. The heat P taken from the 
N electrode is then released in the S reservoir, 4(0) = —Q = 
P+IV, thus the full heat production in both reservoirs is equal 
to the Joule heating, 4(0) +4(L) = IV. 

Now we discuss the boundary conditions. At jc = L, we as- 
sume all functions to be continuous, neglecting a spreading re- 
sistance of the transparent NN' interface: in the diffusive limit, 
this resistance is always small compared to the resistance of 
the N' wire.^^ At the tunnel barrier, x — 0, the function G and 
the matrix current J at the normal (N') and the superconduct- 
ing (S) sides of the junction are connected via the generalized 
boundary condition due to Nazarov,^"* 



1 



rp{r)dr[Gs,GN'] 



2gNRTJo l + ^{{Gs,GN'}-2y 



(12) 



where Rj is the barrier resistance and p(r) is the distribu- 
tion of the transparencies of the conducting channels of the 
barrier (Jq rp(r)iir =1). Assuming the absence of highly 
transparent channels with F ^ 1 and considering p (F) to be 
locaUzed around a small value of F ^ 1 (tunnel limit), we can 
neglect the anti-commutator term in Eq. (fTSI) . thus arriving at 
the Kupriyanov-Lukichev boundary condition^^ alx = 0. 



JN' = i28NRT)-'[GsA 



The boundary conditions for the functions g and G'^ at the 
tunnel barrier follow from Eq. ( fT3] l, 

JN'^iW/^o)[gs,gN'i (14a) 
J§, = {W/^o)[Gs,GN'f. (14b) 

Here the tunneling parameter W is defined as 

W = R{^q)I2Rt = (3^o/4^)r » F, (15) 

where /?(<^o) = ^qEn^ is the resistance of the N' lead per length 
^0- It has been shown in Refs. '26' and '27 that it is this quan- 
tity, rather than the barrier transparency F, that plays the role 
of the transparency parameter for diffusive tunnel junctions. 
Below we consider the case W ^ 1, which corresponds to the 
conventional tunneling limit. 

The N and S electrodes are assumed to be equilibrium reser- 
voirs with unperturbed spectral characteristics and equilib- 
rium quasiparticle distributions. 



On ~ 0, f±N - 



9s = arctanh - 



1 / E + eV E-eV 
- tanh — ± 



2Tn 



/- 



+s 



2Ts 



: tanh - 



f- 



2Tn 
= 0, 



(16) 
(17) 



where and are the temperatures of the N and S reser- 
voirs, respectively. 

Using the parametrization in Eq. (|5]l, we rewrite Eq. ^ as 
the Usadel equation^^ for the spectral angle 9{E,x), 



iT)df:e =2Esmhe. 



(18) 



Here and below we omit the subscript N' for the functions f± 
and in the N' lead. The boundary conditions for Eq. (fTSb 
follow from Eqs. ( fTSI l and ( I14ab . 



a,0|^^o = (2H^/^o)sinh(0o-0s), 



(19a) 
(19b) 



where Oq denotes the value of at = 0. 

The kinetic equations for the functions f± follow from 
Eq. (HI and have the form of conservation laws for the spectral 
currents I±, 



D±dxf± ^ I±~ const. 



(20) 



N'\ 



(13) 



The continuity of the distribution functions at the N'N inter- 
face implies the conditions f± {E,L) — f±N{E). The boundary 
conditions at the SN' interface follow from Eqs. ( I14bb .^t- 

gw/_(£)=G^(£)/_o(£), (21a) 
gNl+iE) = G+{E)[UoiE)-f+siE)], (21b) 

where the subscript denotes the function values at x = and 

G±(£) = Rt\NsNn' TM±M± ), (22) 
A^(£') =Re(cosh0), M+(£') +/M"(£') = sinh0. 

The function A^(£') is the density of states (DOS) normalized 
to its value Np in the normal state; the quantities G± can be 
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interpreted as spectral conductances of the tunnel barrier for 
the probability (+) and electric (-) currents, respectively. At 
large energies, \E\ ^ A, when A^(£') approaches unity and the 
condensate spectral functions (E) turn to zero at both sides 
of the interface, the conductances Gj{E) coincide with the 
normal barrier conductance, R^^. Within the subgap region 
\E\ < A, Gj{E) turn to zero, which reflects blocking of the 
probability current due to full Andreev reflection. 

In the superconducting reservoir, the density of states 
Ns{E) and the condensate spectral functions M'^{E) read. 



NsiE) 
M,{E) 



|£|Q(|£|-A) 

V£2 - A2 ' 
A0(A-|£|) 



MUE) 



(23a) 



-, (23b) 



VA2-£2 ^ ^ ' V£2-A2 

where @{x) is the Heaviside step function. 

B. Solution of the Usadel and kinetic equations 

Generally, the solution of the Usadel equation for a N' lead 
of finite length can be found only numerically. However, in 
the case of a low-transparent tunnel barrier, W <C 1, the spec- 
tral angle is small, <C 1, for all essential energies, which 
enables us to linearize Eqs. (fTsT i and ( |19bl l. 



2Ee, 

= (2W/<^o)(0ocosh0s-sinh0s). 
The analytical solution of these linearized equations 



(24a) 
(24b) 



0(£,x) = 0o(£) 
0o(£) 



sinh[^Af(L-x)/^o] 
2Wsinh05 



kN ■ 



feiV coth (A:a,L/ (^0 ) + 2 W cosh 0s ' 



(26) 



was found to differ from the numerical solution of the exact, 
nonlinearized Usadel equation by less than 1 % for reasonable 
values of W < 10^^. Note that in our approximation we keep 
a small term of the order of W in the denominator of Eq. ( |26] l 
which prevents divergence of Oq at the gap edge, E = A, and 
thus provides a good agreement with the numerical solution 
in the vicinity of this "dangerous" point. 

The analytic solution of the kinetic equations ( |20l ) with cor- 
responding boundary conditions ( 1211 1 is 



f-^f-N- 
f+ = f+N - 



Rj{E)+RNa-{Qy 
{f+N - f+s)RNa+{x) 
R+{E)+RNa+{Q) ' 



(27a) 
(27b) 



a±(x)= / —Di\E,x'), 



where Rj (E) — [Gj {E)\ are spectral resistances of the tun- 
nel barrier,2^ and R^ — Lg]^^ is the normal resistance of the N' 
lead. In Eqs. ( |27] |. we used the relation 



RN/RT=2WiL/^o), 



(28) 



following from the definition of the parameter W in Eq. ( fTSl ). 
In a typical experimental situation, the tunnel resistance dom- 
inates, Rt ^ Rn, therefore the functions f± are always close 
to the equiUbrium distributions f±N in the N reservoir 



C. Electric current 

The electric current is given by the equation obtained by 
combining Eqs. (|6]l, (I21al ) and ( I27al ). 



1 

/= - 



f-NiE) 



eJo R^{E)+R^{E) 



dE, R^{E)=RNa-{0). (29) 



A similar result has been obtained for a NINIS structure in 
Refs. I29I and |2l]; it differs from Eq. ( |29] ) by an additional 
tunnel resistance of the NIN interface in the denominator 
In the spirit of circuit theories for mesoscopic superconduct- 
ing structures, 2^'24 jj^jg equation can be interpreted as "Ohm's 
law" for the spectral current induced by the effective potential 
/_Af in the series of the tunnel resistance Rj and the resistance 
R]^ of the N' lead renormalized by the proximity effect. 

The current in Eq. ( |29l ) involves contributions of both the 
single-particle and the two-particle (Andreev) currents. It is 
useful to discuss these two components separately. To this 
end, we divide the total range of energy integration into two 
regions, £ > A and < £ < A, and take into account that the 
superconducting DOS A^s(£') = at < £ < A and the spectral 
function (£) = at £ > A, 



I = h+h 



f-NiE) 



eJA Rt{NsNn>)-'+R^{E) 



f-N[E) 



eJo Rt{MsM^,)-^+RJ,{E) 



dE 



dE. (30) 



The main contribution to the current Ii comes from the pro- 
cesses of single-particle tunneling. Besides, Ii contains small 
proximity corrections due to deviations of the DOS A^^?/ and of 
the diffusion coefficient D in the N' lead from their unper- 
turbed values Nn = Dn = 1 . Physically, these deviations are 
due to the partial Andreev reflection at the energies above the 
superconducting gap and therefore rapidly decay as the energy 
increases. Neglecting this small effect, we obtain the formula 



1 



f-N{E) 



e J A RtN. 



-1 



Rn 



-dE, 



(31) 



which describes the single-particle current in the NN'IS struc- 
ture. At the subgap voltages, eV < A, this current tends to 
zero exponentially at small temperatures, Tn <C A. At large 
voltage, eV ^ A, the current /j approaches an Ohmic depen- 
dence with the deficit current arising from the contribution of 
the N' lead to the net junction resistance R = Rj +Rn, 



J y , , rA V2 Rti 
R eR r Rj 



eR 



(32) 



Finally, neglecting the small contribution Rn to the junction 
resistance and rewriting f^N in terms of the Fermi function of 
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the N reservoir, n^iE) = [1 +exp(£'/7)v)] ^ we arrive at the 
standard formula of the tunnel theory ,22 

h = -^[ Ns{E){nN{E-eV)^nN{E)]dE. (33) 

Within the same approximations, the Andreev current is re- 
duced to the following form. 



Ia = 



8±iE) 



2WA^ 



f-N{E) 



sRt Jo sl + ig- 

sinhjS ±sinj3 
coshjS — cosjS 



-2WEy VA2-£2 



£(A2-£2) 
2A ' 




At large voltage, eV ^ A, the Andreev current approaches 
a constant value I^xc ~ (\/2WA/e/?)ln(V2/W) (excess cur- 
rent); for long junctions, L ^ ^q, it is much smaller than the 
single-particle deficit current in Eq. ( [32] i. Thus in this limit the 
net electric current, /i +Ia, always exhibits a deficit current. 
The Andreev current in NIS structures was first calculated mi- 
croscopically by Hekking and Nazarov'^ (see also Ref. I3 ll) 
and VoUcov et alM Note that in our consideration we neglect 
possible pair-breaking factors (like magnetic impurities) and 
damping of quasiparticles in the S region due to inelastic in- 
teractions. For this reason, our results concerning Andreev 
current may differ from that by Volkov et alr^^ especially at 
small eV comparable with corresponding relaxation rates. 

We would like to notice that the Andreev current does not 
depend on the N' lead length L as long as L ^ ^o- In this case, 
the magnitude of the Andreev current at eV ~ A can be esti- 
mated from Eq. dUi as Ia ^ WA/eRr = AR{^o)/2eR^. This 
reproduces the result of Hekking and Nazarov'^ and Volkov 
et alM As energy decreases, the spectral density of the An- 
dreev current [the integrand in Eq. (|34] |1 diverges as £'^'/2 
until E reaches the small Thouless energy En ^Ti/L^, which 
plays the role of a cut-off factor Such behavior of the An- 
dreev current was first discovered in Ref. 16 using the di- 
agrammatic methods in the tunnel Hamiltonian formalism. 
In the limit of a short junction, L ^ ^q, when the prox- 
imity effect and the Andreev current are suppressed by the 
N reservoir, we recover the result of Bardas and Averin,'^ 
Ia - AR{L) /leRj. = M^N/leR^. 



D. Energy current 

The energy current can be obtained upon combining 
Eqs. dU, (l2lf] i. and (l27b] i, 

e^Jo R+[E)+R+{E) ' ™ +v ; 



(35) 



This expression is quite similar to Eq. ( |29l l for the electric 
current and has the same physical interpretation: the spectral 
probability current flowing through the series of the tunnel 
and normal resistances is determined by Ohm's law for the 
effective potential difference f+N ~ f+s- 



First, we note that the energy integration in Eq. ( l35T l is actu- 
ally confined to the interval E > A since the conductivity Gj 
turns to zero (and, correspondingly, Rj ^ °o) at Q < E < A. 
Thus the Andreev energy current Qa is identically zero; phys- 
ically, this corresponds to the fact that the quasiparticle prob- 
abiUty current /+ is completely blocked in the subgap energy 
region due to full Andreev reflection. 

Neglecting the proximity corrections to the spectral func- 
tions, i.e., assuming Nj^i ~ D+ = 1 and Mi = 0, we obtain a 



dE. (34) simplified form of the single-particle energy current. 



Qi 



^Ei±im^MEldE. 



RtN^ 



Ra 



(36) 



Finally, omitting the contribution Rn of the normal lead to 
the total resistance and expressing /+ in terms of the Fermi 
functions we arrive at the standard form for the energy current, 

Qi=-^r Ns{E)E[nN{E-eV)-ns{E)]dE, (37) 



?^R 



T J- 



where ns{E) = [1 +exp{E/Ts) 
S reservoir. 



is the Fermi function of the 



E. Heat current 

The heat current out of the normal metal reservoir (cooling 
power) can now be obtained from the above expressions for 
the electric and the energy currents, Eqs. ( [30l l and ( [35] l, using 
Eq. (fTTT l. As follows from Eq. (fTTl i. the Andreev heat current 
to the normal reservoir is nonzero giving a negative contribu- 
tion Pa to the cooling power. 



P = Pi+Pa, Pa = -IaV. 



(38) 



From this equation we see that the heat current out of the nor- 
mal metal is affected by the Joule heating generated by the 
Andreev current Ia- This is due to the fact that the Andreev 
current is fully dissipated in the normal metal. 

Using the tunnel model formula ( [33] ) for the electric cur- 
rent and Eq. ( l37T i for the energy current, we arrive at the well- 
known form for the cooUng power^^ 

Pi = ^ r Ns{E){E-eV)[nN{E-eV)-ns{E)]dE. 

(39) 

This equation is widely used when fitting the experimental 
data on electron cooling. Such an approach is valid as long as 
the Andreev contribution to the electric current is negligibly 
small, i.e., at moderately high temperatures. As noted above, 
at low temperatures, the single-particle processes are expo- 
nentially suppressed in the subgap voltage region, where the 
effect of Andreev current on electron cooling becomes essen- 
tial and must be taken into account. 



III. PLANAR NN'IS MODEL 

In this section we present an extension of the approach 
developed above to the more realistic case of a sandwich- 
type tunnel junction with a thin-film N' lead as sketched in 
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Fig-HIb). This situation is more complex; however, it is possi- 
ble to reduce this problem to the ID case by formulating effec- 
tive boundary conditions at the junction, following a method 
suggested by Volkov^^ and Kupriyanov^^ 

In the general three-dimensional case, the Keldysh-Usadel 
equations (Ull and the boundary condition Eq. ( fTsT i read 



[a,£,G] =/'DVJ, J^GVG, 



(40a) 
(40b) 



where n is a vector normal to the insulator layer. In Eq. ( I40bb 
all functions are taken at the sides of the barrier. 

We suppose the size of the planar junction Lh to exceed the 
coherence length, Li, » ^q, and the thickness of the N' lead to 
be much smaller than the coherence length, d ■^£,0. Then the 
function G in the left-hand side of Eq. ( I40al i is approximately 
constant within the normal metal bank above the junction . 
Upon integration of this equation over the volume of the nor- 
mal metal bank, transforming the volume integral in the right- 
hand side into a surface integral, and using the boundary con- 
dition Eq. ( |40b| ) at the tunnel barrier, we obtain the effective 
boundary condition for the ID Keldysh-Usadel equations in 
the N' lead, 

Sbd[a,E,GN,] = iV{SJo- Sh{W /^o)[Gs,Gn']}, (41) 

where all functions are taken at the sides of the barrier In 
Eq. dTTT i. S is the cross-section area of the N' lead, d is the 
lead thickness, Sh is the area of the junction [see Fig. [Ub)], 
and Jq = GffidxGf^i is the value of the matrix current in the N' 
lead at the cross-section adjoining the junction (i.e., atx = 0). 
A similar result has been obtained in Ref. [ss] for the case of a 
planar SIS junction. Equation W\\ can be rewritten as 

[g.E.Gn'] = 2/A{(^oV'Li)/o- W[Gs,Giv']}, (42) 

where 



W^W{^Q/d)^[3^^,/AU)T 



(43) 



is the effective tunneling parameter. Note that for thin-film 
planar junctions this parameter is much larger than the ID tun- 
neling parameter W by the ratio E,Q/d^ 1 . 

As long as E,q <C if, and E,qJq ~ W , the first term in the right- 
hand side of Eq. (|42] | can be assumed to be the smallest one 
and thus neglected. However, this is only true for the Green 
component of Eq. (|42]) . 



(44) 



whereas for the Keldysh component the diagonal part of the 
left-hand side of Eq. (l42l l turns to zero, and therefore the 
boundary condition for the diagonal part of reads 

-/^ = (lV//^o)[Gs,Gw']^, (45a) 
Wf ^ W{Lh/d) = W{Sh/S) = W{Lh/^Q) » W. (45b) 

The enhancement of the parameter Wf with respect to W re- 
flects decrease of the tunnel resistance Rt compared to its 



value in the ID case, due to increase of the junction area from 
S for the ID geometry to Sh for the planar geometry (provided 
the barrier transparency is equal for both cases). 

In terms of the spectral angle 9 in the N' lead, the boundary 
condition (|44] | has the form k]^ sinh 9q — 2W sinh(05 — 9o) and 
can be solved explicitly for the boundary value 9o, 



9q = arctanh 



2Wsmh9s 
kl + 2W cosh 9s' 



(46) 



Equation ( l46b results in a DOS minigap in the normal bank of 
the junction. To first order in W this minigap is equal to 



' 2WA < A, 



(47) 



and the spectral angle 9o is given by the BCS-like formula 
00 ~ arctanh (Ag/ii) at small energies, E ^ Ag. The spatial 



dependence of the spectral angle in the N' lead obeys Eq. 
with 00 defined in Eq. ( l46b . 

We note that a similar result was found for short SINIS 



junctions 



36,37,38 



The analogy between the NIS sandwich and 



a short SlNlS junction can be clearly seen from the mapping 
method, similar to the one used in electrostatic problems. In- 
deed, at the top surface of the N bank the boundary condition 
reads 39 /dn = 0. To ensure this condition, we add a mirror 
image of the NIS sandwich to the top surface of the N layer, 
thus arriving to the problem of a SlNlS junction with a normal 
metal interlayer of thickness 2d. 

The boundary condition ( 145 a| ) for the distribution functions 
is similar to Eq. ( I14bb in the ID case, with the substitution 
W ^Wf. As follows from the definition of the parameter Wf 
in Eq. ( |45b| i, the ratio of the normal and tunnel resistances is 
similar to Eq. ( |28] l, Rn/Rt = 2WfL/^o; therefore, the distri- 
bution functions in the planar geometry, being expressed in 
terms of the spectral resistances, coincide with the result for 
the ID case, Eqs. dZTl i. As a result, equations ( l30l l and ( l35l l for 
the electric and energy currents hold their form for the planar 
geometry, however with different tunnel resistances Rj. 

We note that within the main approximation in W, the 
spectral density of the Andreev current is nonzero only in- 
side the minigap, E < Ag. In this energy region, the spec- 
tral functions and M^, are approximately equal to — 1 and 

— Ag/(Ag — respectively. Using Eq. dSOl l and neglect- 

ing small contribution of the N' lead to the net resistance, we 
obtain a simple expression for the Andreev current at eV ^ Ag 
in the planar NIS junction. 



1 eV f^s 

Ia^— tanh — / M M ,dE 
eRr 2Tn Jo 



TzA,, , eV 
— tanh 

2Tn 



2eRi 



(48) 



IV. RESULTS AND DISCUSSION 

In our numerical calculations and analysis, we use exact ex- 
pressions for the electric, energy, and heat currents, Eqs. dSOb . 
(|35] l, and ( fTTT i. taking into account both the proximity correc- 
tions and the resistance of the N' wire. Although these effects 
give small contributions to the energy and electric currents 
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eWA(r) eVI^(T) 

FIG. 2; (Color online) Cooling power versus bias voltage at IV = 0.5 • 
10"^ (a) and IV = 2 • 10"^ (b) for different temperatures: T = Q.lTc 
(red line), T = Q3Tc (blue line), and T = Q.STc (black line). Solid 
lines represent full cooling power, dashed lines were computed at 
/a = 0. 

separately, the cooling power P, being a relatively small dif- 
ference of the energy current and the Joule heat, is very sen- 
sitive to small details of the charge and energy transport. We 
will focus on the case of a planar junction which is the most 
adequate model of a real experimental setup. In what follows, 
we assume the temperatures of the N and S reservoirs to be 
equal, Tn = Ts = T. 

We note that the tunneling parameters W and W, accord- 
ing to their definition in Eqs. (fT5l),(|43]l, are temperature- 
dependent, since the coherence length £,q increases with tem- 
perature as A^'/^(r). This variation is important at high 
enough temperatures and it was taken into account in our cal- 
culation scheme, although at low temperatures, T <C 71, when 
the cooling effect becomes apparent, one may neglect the tem- 
perature dependence of A. In what follows, we assume the 
quantity W and its effective value W in planar junctions to 
be taken at T = 0, allowing for their temperature dependence 
in Eq. (|46] l for the spectral angle by means of corresponding 
temperature-dependent factors. In order to keep a common 
scale of the cooling power P calculated for different T and W, 
we normalize P to the ratio A^{0) /e^Rro, where Rto is the 
junction resistance at a fixed value W = 10^^ of the tunneling 
parameter. Relying on typical sizes of the experimental sam- 
ples, we assume L — Lj, — lO^o. where the coherence length 
^0 is taken at T = (for Al-based film structures, its value is 
about 100 nm). 

Now we proceed to the discussion of our results. The effect 
of the Joule heat generated by the Andreev current 7^ on the 
cooling power is illustrated by voltage dependencies Piy) in 
Fig. |2] where the solid curves depict full cooling power, and 
the dashed curves were calculated at 7^ = 0. For a highly- 
resistive tunnel junction \W — 0.5 • 10^^, Fig.|2ja)], the heat- 
ing effect due to Andreev current is negligibly small. For 
smaller junction resistance \W — 2- 10^^, Fig.|2b)], the heat- 
ing effect essentially modifies the result; in particular, at low 
enough temperature, T = 0.1 r^, it makes P{V) negative at all 
voltages. This is due to the fact that for phase-coherent dif- 
fusive proximity systems, the two-particle contribution to the 
subgap transport is anomalously strong at low energies. 

As it is obvious from Fig.|2l the cooling power approaches 
a maximum at a certain optimal bias voltage Vopi which de- 
pends on both the temperature and the tunneling parameter. 




0.0 0.5 1.0 0.0 0.5 1.0 



r/r T/T 

c c 

FIG. 3: (Color online) Temperature dependencies of the optimum 
bias Vopi{T) (a) and of the cooling power P(T) at optimum bias (b), 
for different values of the tunneling parameter: W = 10^"* (black), 
W = 10-^ (blue), and W = lO^^ (^ed). 




FIG. 4: (Color online) Dependencies of the optimum bias V„pr{W) (a) 
and of the cooling power P{W) at optimum bias (b) on the tunneling 
parameter, for different temperatures: T = 0.1 (red), T = 03Tc 
(blue), and T = O.STc (black). 

It is interesting to note that the dependence Vopt{T) is almost 
universal within a wide range of the tunneling parameter, as 
shown in Fig. lUa). AtT > Q.75Tc, the optimal bias formally 
turns to zero which means that at these temperatures the cool- 
ing power becomes negative for all voltages. Existence of the 
upper limiting temperature for the cooling effect is explained 
by the increase in the number of thermally excited quasiparti- 
cles which produce enhanced Joule heat. As the temperature 
decreases, the optimal bias rapidly increases and approaches a 
value rather close to the energy gap A(r). Simultaneously, the 
cooling power at optimal bias, P[Vopt{T)], first increases and 
approaches a maximum at T « (0.4 0.5)rc-, see Fig. (Sjb). 
Then, at lower temperatures, the Joule heat due to Andreev 
processes causes the cooling power to decrease. At a certain 
temperature 7;,„„, the cooling power tends to zero, which de- 
fines the lower limiting temperature for the cooling regime. 
As follows from Fig.[3jb), the temperature T,„i„ increases with 
the tunnel parameter, approaching 0.24-Tc for W = 10^^; this 
is because the Andreev cuiTent and the associated Joule heat 
increase with the junction transparency more rapidly than the 
single-particle cooling power. At temperatures just above 
T„i„, the function P{V) at small applied voltage is negative; 
thus, the electron temperature is expected to increase first with 
the bias due to Andreev current heating before it decreases 
due to the single-particle cooling effect. This phenomenon 
has been observed in experiments^^ at very low temperatures. 
The dependencies of the optimal bias on the tunneling pa- 
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rameter W are plotted in Fig. Ufa). In accordance with the 
above-mentioned universaHty of the curves Vopt{T) for dif- 
ferent values of W, the dependence Vopt{W) is rather weak 
at W < 10^2 Within this region the cooling power at opti- 
mal bias, P[Vop,(W)], linearly increases with W, as shown in 
Fig. SJb), which is expected when single electron tunneling 
dominates. For larger values of the tunnel parameter, the An- 
dreev current heating dominates over the single-particle cool- 
ing and leads to a rapid decrease of the cooling power, which 
tends to zero at a certain onset point, as seen from Fig.HJb). 
As the temperature decreases, the role of Andreev processes 
becomes more important, therefore the onset shifts towards 
smaller values of W. 



V. EXTENSION TO SINIS JUNCTION 

As noted in the Introduction, in most experiments the re- 
frigerator is arranged as a double -barrier SINIS junction,'^ 
where the S electrodes are massive reservoks, and N is a nor- 
mal metal strip. Generally, in such structures the charge and 
energy transport is due to multiple Andreev reflections (MAR) 
of quasiparticles from the NS interfaces.^^'^^ During every 
passage across the junction, the electrons and the retro-reflec- 
ted holes gain an energy eV, which allows them eventually to 
overcome the energy gap and to escape into the S reservoirs. 
This results in a strong quasiparticle nonequilibrium charac- 
terized by intense electron heating within the subgap energy 
region, which has been detected in the experiments.'*" From 
this point of view the cooling effect observed in SINIS junc- 
tion looks at a first glance somewhat surprising. 

However, the inelastic scattering processes impose strong 
limitations for the existence of the MAR regime: in order to 
provide quasiparticle diffusion through the whole MAR stair- 
case, from —A to A, the quasiparticle dwell time in the N lead, 
Tf/, must be smaller than the inelastic relaxation time Tg (for 
details, see Ref. 21). In typical cooling experiments on SI- 
NIS junctions with low-transparent SN interfaces,'^ the dwell 
time greatly exceeds t^, which prevents accumulation of the 
quasiparticle energy gains and thus destroys the MAR regime. 
Correspondingly, the distribution functions in the N lead be- 
come close to local-equilibrium ones. 



f±iE,x) 



1 



tanh- 



-eV(x) , E-eV(x) 
' -itanh ^ ' 



2Tn 



2Tn 



(49) 



where V{x) is the voltage at the given point x. The variations 
in V{x) are of the order of V{Rn/R), i.e., negligibly small 
compared to the applied voltage V which mainly drops at the 
tunnel barriers. This implies that the distribution functions are 
close to equilibrium functions in a normal reservoir In this 
case the SINIS junction behaves as two NIS junctions con- 
nected in series through the equilibrium normal reservoir. As 
for the spectral angle, in long junctions, L ^ it can be ap- 
proximated by the solution of the Usadel equation for a semi- 
infinite NIS structure;'^ simultaneously, this solution is also 
a good approximation to our solution for a long NN'IS struc- 
ture. 



From this we conclude that our results can be applied to 
the description of electron cooling in SINIS structures with 
large quasiparticle dwell times. Similar modelling of a SINIS 
junction by a series of two NIS junctions has been used in 
Ref.i32ifor the calculation of the differential conductance. 



VI. SUMMARY 

We have developed a quantitative theory of charge and heat 
transport in one-dimensional and planar NN'IS tunnel junc- 
tions and studied the effect of electron cooling in such struc- 
tures. We extend the microscopic approach by Bardas and 
Averin,'^ originally applied to constriction-type junctions, to 
structures of arbitrary length and thin-film geometry used in 
practice for microcooler fabrication. We found that the contri- 
bution of two-particle (Andreev) current to the Joule heat gen- 
erated in the normal reservoir noticeably modifies the cooling 
effect, especially at low temperatures and/or in rather trans- 
parent junctions. The interplay between the Andreev cur- 
rent heating and the single-particle cooling, whose intensity 
rapidly decreases with temperature, determines the lower lim- 
iting temperature T,„i„ for the cooling regime. When the trans- 
parency of the NIS interface increases, the Andreev processes 
play a more essential role, therefore the temperature r„„„ in- 
creases. At high temperatures, the cooling regime is confined 
by the enhancement of the Joule heat due to thermally ex- 
cited quasiparticles; the maximal cooling temperature is about 
0J5Tc, being almost independent of the junction resistance. 
As a result, the cooling effect persists within a specific tem- 
perature interval and approaches a maximum at the tempera- 
tures (0.4^0.5%. 

We pay special attention to the analysis of the optimum bias 
voltage V(,pt{T,W), at which the cooling power approaches a 
maximum for given temperature T and the tunneling param- 
eter W. We found that Vopt exhibits a virtually universal tem- 
perature dependence for different values of the tunneling pa- 
rameter and approaches values close to the energy gap as long 
as the temperature decreases. The cooling power at optimum 
bias voltage first increases linearly with W until the Andreev 
current heating abruptly suppresses the cooling regime. 

We discussed the applicability of our results to the descrip- 
tion of the cooling effect in SINIS junctions. We show that 
such a double-barrier structure can be modeled by a series 
of two independent NIS junctions, provided the quasiparti- 
cle dwell time inside the junction greatly exceeds the inelastic 
relaxation time. This condition, which is usually satisfied in 
cooling experiments,''* enables one to extend the theory pre- 
sented here to the case of the SINIS microcoolers. 

From our considerations we conclude that the Andreev cur- 
rent is one of the most serious factors of limitation of the elec- 
tron cooling efficiency. In order to reduce this factor, one 
should address materials in which the proximity effect and, 
correspondingly, the Andreev current are strongly suppressed. 
A first guess to such materials can be ferromagnets.'" A quan- 
titative analysis of the cooling effect in FIS junctions will be 
presented elsewhere. 
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